Theory of prospective tetrahedral perovskite ferroelectrics 
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Using first-principles methods, we predict the energy landscape and ferroelectric states of double 
perovskites of the form AA'BB'Oe in which the atoms on both the A and B sites are arranged in rock- 
salt order. While we are not aware of compounds that occur naturally in this structure, we argue 
that they might be realizable by directed synthesis. The high-symmetry structure formed by this 
arrangement belongs to the tetrahedral Fi3m space group. If a ferroelectric instability occurs, the 
energy landscape will tend to have minima with the polarization along tetrahedral directions, leading 
to a rhombohedral phase, or along Cartesian directions, leading to an orthorhombic phase. We find 
that the latter scenario applies to CaBaTiZrOe and KCaZrNbOe, which are weakly ferroelectric, 
and the former one applies to PbSnTiZrOe, which is strongly ferroelectric. The results are modeled 
with a fourth- or fifth-order Landau-Devonshire expansion, providing good agreement with the 
first-principles calculations. Computations of zone-center soft modes are also carried out in order to 
characterize the polar and octahedral-rotation instabilities in more detail. Prospects for synthesis 
of ferroelectric materials belonging to this class are discussed. 

PACS numbers; 77.80.bg, 77.84.Bw, 81.05.Zx, 71.15.Nc 



I. INTRODUCTION 

Recent interest in novel materials has been stimulated 
by unprecedented advances both in experimental materi- 
als synthesis and in first-principles computational meth- 
ods for predicting materials properties. Attention has 
focused in particular on functional materials that can 
be driven between various structural or electronic phases 
having distinct properties, as for example by the appli- 
cation of electric fields, magnetic fields, or strain. 

Ferroelectric perovskites constitute a subset of these 
interesting compounds, with their switching behavior 
providing potential applications in non- volatile memories 
and their piezoelectric properties making them attractive 
as actuators and sensors. The perovskites are also of con- 
siderable interest in the search for multiferroics having 
strongly coupled polar and magnetic properties. Recent 
work has shown that many perovskite properties can be 
tuned through the application of epitaxial straiuii"— 

Ferroelectric perovskite oxides can generally be clas- 
sified into those derived from cubic symmetry, meaning 
that the actual or putative high-temperature symmetric 
phase is cubic, or those derived from tetragonal symme- 
try, meaning that the high-symmetry phase is tetrago- 
nal. Most well-known ferroelectrics, including BaTiOs, 
KNbOa, PbTiOs, and BiFeOa belong to the first class, 
while some layered ferroelectrics, such as SrBi2Ta209,4"- 
belong to the latter. It is possible for the compositional 
ordering in tetragonal layered systems to break inversion 
symmetry, as for example in "tricolor" super lattices i^""— 
To our knowledge, however, there are no perovskite ox- 
ides in which the chemical composition corresponds to 
a tetrahedral high-symmetry structure. While some bo- 
racites (M3B7O13X, where M is usually a divalent metal 
and X is usually a halogen) such as Ni3B70i3l are realiza- 
tions of this kind of tetrahedral ferroelectric system^""— 
it would be very interesting to see the same symmetry 



class represented in the better-known cases of perovskite 
oxides. 



The present work is motivated by the idea that per- 
ovskite oxides having a tetrahedral high-symmetry struc- 
ture might be realized experimentally and might have in- 
teresting ferroelectric or other physical properties. With 
this in mind, we have carried out a computational study, 
based on first-principles density-functional calculations, 
of perovskites having a tetrahedral compositional sym- 
metry. In particular, we focus on AA'BB'Og double per- 
ovskites in which both A and B sites exhibit rock-salt 
order. While rock-salt order on the B site of a double 
perovskite is common, it is quite rare on the A site, where 
a 50% mixing of two atoms typically leads to (001) lay- 
ered ordering if it orders at all. Nevertheless, even if the 
tetrahedral symmetry is not realized in the equilibrium 
phase diagram, we shall argue that experimental routes 
to the directed synthesis of such double-rock-salt tetra- 
hedral AA'BB'Og perovskites may be available. 

The paper is organized as follows. Sec. |TT] introduces 
the structure of the AA'BB'Og perovskite materials, de- 
scribes the possible symmetry-determined directions of 
the polarization, and discusses possible domain types 
and their symmetries. We then briefly detail our the- 
oretical methods in Sec. IIIIl Sec. |IV] describes the re- 
sults of our first-principles calculations and presents a 
comparison with a simple empirical model. We also dis- 
cuss the effects of rotations of oxygen octahedra, and 
present some results concerning a Mn-containing double 
perovskite in which magnetic ordering is also an issue. 
Finally, in Sec. |V] we summarize our work and present 
our conclusions. 



2 




FIG. 1: (Color online) Structure of AA'BB'Oe double per- 
ovskites, with A (red circles) and A' (magenta squares) atoms 
forming one rock-salt framework, and B (blue diamonds) and 
B' (green triangles) atoms forming a second interpenetrat- 
ing one. Oxygen atoms are not shown, although one oxygen 
octahedron is outlined to clarify their role in the structure. 
Dashed lines illustrate the tetrahedral point symmetry. 



II. TETRAHEDRAL AA'BB'Oe DOUBLE 
PEROVSKITES WITH ROCK-SALT ORDER 



Cubic Pm3m parent structure 
(a). 




(lll)^i?3m 
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Tetrahedral F43m parent structure 
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FIG. 2: (Color online) Symmetry-determined possibilities for 
directions of energy minima in the space of polarization, rep- 
resented by red dots. Fig. (a)-(c) cubic composition leading 
to high-symmetry Pm3m structure, and (d)-(f) tetrahedral 
composition leading to high-symmetry F43m structure. 



We consider here a class of perovskite ferroelectrics 
whose compositional symmetry, and thus the high- 
temperature symmetric phase, is tetrahedral, instead of 
cubic or tetragonal. The simplest way to arrive at tetra- 
hedral symmetry in the perovskite system is to populate 
both A and B sites with two different kinds of atoms (A 
and A', and B and B', respectively) arranged in rock- 
salt (three-dimensional checkerboard) order. The crys- 
tal chemistry of the perovskites, while preferring rock- 
salt ordering for B sites, resists the same on the A sites, 
where layered (001) ordering is preferred insteadJ^iii To 
our knowledge, no perovskite oxides exhibiting simulta- 
neous A- and B-site rock-salt ordering has been reported. 
Nevertheless, we shall investigate their properties theo- 
retically here. 

The structure in question is illustrated in Fig. [TJ in 
which oxygen atoms are suppressed for clarity. It can be 
seen that the point symmetry of each atom is tetrahe- 
dral, and since there is only one formula unit per primi- 
tive cell, this also establishes the crystal point group as 
tetrahedral, with the F43rn space group (216). This then 
identifies the parent high-symmetry structure, and can be 
compared with the Pm3m high-symmetry structure that 
characterizes most ordinary perovskite ferroelectrics. 

In such a tetrahedral ferroelectric, the symmetry of 
the energy landscape E{P) is such that E will be sta- 
tionary with respect to the direction of P for P along 
(111), (ill), and (100) directions. If there are local min- 
ima along these directions, they correspond to the rhom- 
bohedral R3m space group (160) in the first two cases, 
or the orthorhombic Imm2 space group (44) in the last 
case. This is illustrated in Fig. [21 where panels (a-c) il- 



lustrate the usual case of cubic perovskites, while panels 
(d-f) show the corresponding possibilities in the tetrahe- 
dral perovskite system. If the rhombohedral directions 
are energetically favored, as in panels (d-e), then there 
will be two distinct types of rhombohedral domains with 
different energies, and the dependence of the total en- 
ergy on P for P = Pn along the body diagonal n will 
display an asymmetric double-well potential. This is in 
sharp contrast with systems like A2BB'06 or AA'B206, 
in which only one sublattice has rock-salt order; the sym- 
metry is then that of panels (a-c), and the magnitude of 
the polarization does not change when the polarization 
is reversed during a domain switching event. 

Even if these perovskites with double rock-salt do not 
exist as equilibrium phases, it may be possible to synthe- 
size them using atomic-layer epitaxy techniques. In order 
to access the desired compositional order, this would need 
to be done by epitaxial growth on the (111) perovskite 
surface orientation. As a proof of principle, Rijnders, 
Blok and Blank have recently reported such a growth 
of SrCaTiMnOg and CaBiTiFeOe films on LaAlOg and 
SrTiOs substrates. Interestingly, a symmetry analysis 
shows that a (111) uniaxial strain, of the kind that results 
from a lattice mismatch in this kind of epitaxial geom- 
etry, would convert the system to a polar space group 
by selecting out one of the body-diagonals as special, re- 
sulting in a polarization even when none is present in the 
high-symmetry tetrahedral structure. A related effect, in 
which improper ferroelectricity can be induced by octa- 
hedral rotations around a body-diagonal direction, will 
be discussed in Sec. IIVDI 
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TABLE I: Compounds for which calculations have been car- 
ried out. Second column gives the acronym that we use to 
identify the compound. Lattice constants are calculated in the 
high-symmetry FiSm structure and are reported in terms of 
an effective 5-atom cubic cell dimension. "Average" refers to 
the average computed lattice constant of the ABO3, A'BOa, 
AB'Oa, and A'B'Oa parent materials. 



Compound 


Alias 


Lattice constant (a.u.) 






Calculated 


Average 


PbSnTiZrOe 


PSTZ 


7.47 


7.47 


KCaZrNbOe 


KCZN 


7.55 


7.55 


CaBaTiZrOe 


CBTZ 


7.54 


7.53 


KSrTiNbOe 


KSTN 


7.35 


7.33 


KBaTiNbOe 


KBTN 


7.47 


7.48 


SrCaTiMnOe 


SCTM 


7.15 


7.20 



III. COMPUTATIONAL DETAILS 



The calculations are performed using density- 
functional theory (DFT) as implemented in the ABINIT 
code packageji2r~ We use Ceperley- Alder— 1^ exchange- 
correlation in the Perdew-Wang^l parameterization, and 
TrouUier-Martin norm-conserving pseudopotentials con- 
structed using the FHI98PP^^ code. A plane-wave cut- 
off of 70 Hartree is applied. The Brillouin zone of the 
10-atom fee cell is sampled by a 6 x 6 x 6 Monkhorst- 
Pack^^ fc-point mesh, while that for the 20-atom unit 
cell is 6 X 6 X 4. (Energy differences between a low- 
symmetry and a high-symmetry structures are always 
computed using an identical cell and k-point mesh.) A 
stress threshold of 2 x 10~^ GPa is used for cell relaxation, 
and forces on ions are converg ed below 2.5x10^^ eV/A. 
The electric polarization is calculated using the Berry- 
phase approachj^l To avoid any potential confusion in 
the branch choice for the polarization of non-orthogonal 
unit cells, we also estimated the polarization using the 
computed Born effective charge tensors and atomic dis- 
placements, finding good agreement. 

In addition to ground-state relaxation calculations, we 
also use density-functional perturbation theory to com- 
pute the frequencies of the zone-center phonon modes, as 
an aid in identifying polar or nonpolar (e.g., octahedral 
rotation) instabilities. To this end, the phonon frequen- 
cies at the r point of the 10-atom cell (corresponding to 
both F and R points of a 5-atom cubic cell) were calcu- 
lated, and the corresponding soft-mode eigenvectors were 
analyzed for any unstable modes having imaginary fre- 
quency. The plane- wave cutoff and other details were the 
same as for the ground-state calculations. 



IV. RESULTS 

A. High-symmetry states 

As mentioned in Sec. HI rock-salt ordering is not a 
common form of ordering on the A site of mixed per- 
ovskites, and we are not aware of any naturally-occurring 
AA'BB'Oe double perovskites that exhibit rock-salt or- 
dering on both A and B sites. Therefore, our first step has 
been to carry out a theoretical search for potential candi- 
date materials of this kind. Even if the double rock-salt 
ordering is not the ground-state equilibrium structure for 
a given material, it could be a candidate for attempts at 
directed experimental synthesis. Thus, we have carried 
out calculations for a variety of compounds, placing each 
in the high-symmetry F43m structure and minimizing 
the energy with respect to the lattice constant. We look 
for materials that are insulating, and select possible can- 
didates showing a range of ionic sizes, or combinations of 
cationic charges. For example, we consider combinations 
where both parent perovskites have +2 and +A cations, 
or those with one parent perovskite having -f 1 and +5 
cations while the other has +2 and -1-4 cations. This 
search led us to focus on the six candidate materials that 
are listed in Table HI The second column of this Table 
gives the alias by which each compound will be denoted 
in the remainder of the paper ("PSTZ" for PbSnTiZrOe, 
etc.). 

The third column of Table |I] presents our calculated 
lattice constants for the AA'BB'Og double perovskites. 
In order to put these in context, we also calculate the cell 
constants of all relevant parent perovskites in their high- 
symmetry Pmim structure. Then, for each of these six 
compounds, we average the calculated lattice constants 
of the four parent materials (even if some were occasion- 
ally metallic) and present the result as the last column 
of Table HI In each case we find that the cell constant of 
the double perovskite is very close to average of the par- 
ents. For PSTZ, for example, we find that the parents 
PbTiOs, PbZrOs, SnTiOs and SnZrOs have ceh con- 
stants of 7.30, 7.70, 7.23 and 7.66 a.u. respectively. The 
corresponding arithmetic average is 7.47 a.u., essentially 
the same as the calculated value for PSTZ. The fact that 
the agreement is so good for all six cases indicates that 
the volumes of the parent simple perovskites essentially 
determine the cell constants of the double perovskites in 
the rock-salt structure. (Note that the results reported 
for SCTM are calculated in the high-symmetry ferromag- 
netic spin state, even though the antiferromagnetic struc- 
ture is lower in energy; we do this to stay in the spirit of 
reporting the high-symmetry behavior here. More real- 
istic spin structures will be considered in Sec. IIVEI ) 

B. Exploration of polar instabilities 

Next we look for ferroelectric instabilities in these sys- 
tems by checking the high-symmetry i^43m structures 
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TABLE II: Frequency and principal character of lowest zone- 
center mode, and same for next-higher mode, in the Fi3m 
structure for AA'BB'Os materials. 



Phonon I Phonon II 



AA'BB' 


uj (cm ■"") 


char. 


cj (cm ■"") 


char. 


PSTZ 


182j 


O only 


160i 


A', 


KCZN 


195i 


O only 


96i 


A', O 


CBTZ 


U7i 


O only 


75i 


A, 


KSTN 


W3i 


O only 






KBTN 


149 


A, A', O 






SCTM 


172i 


O only 







to see if there are any phonon modes with imaginary 
frequency. Because of the high symmetry, all phonon 
modes at F have three-fold degeneracy. After identify- 
ing and discarding the triplet corresponding to the zero 
mode (uniform translation), we report the lowest rele- 
vant mode frequencies in the second and fourth columns 
of Table HIl We also inspect the mode eigenvectors and 
report the character of these modes in the third and fifth 
columns of the Table. 

We find that KBTN does not exhibit any unstable 
modes, suggesting that it is probably stable (or, at least, 
metastable) in the FA3m structure. Turning now to the 
other five materials, we see that there are two sets of un- 
stable soft modes for the first three materials, while there 
is only one for KSTN and SCTM. Moreover, the most un- 
stable mode always has character only on oxygen atoms, 
indicating that it corresponds to a pattern of octahedral 
tilts or rotations. Since this occurs at the F point of the 
10-atom cell, it corresponds to an i?-point instability of 
the 5-atom parent perovskite. (We shall consider rota- 
tional instabilities further in Sec. IIVDI ) However, if we 
look at the other set of unstable modes for each of the first 
three materials, we observe a large contribution coming 
from the smaller of the A ions, signaling that they are 
polar (i.e., infrared-active) distortions. We henceforth fo- 
cus our attention on an in-depth study of the first three 
materials, namely CBTZ, KCZN and PSTZ. 

As discussed in Sec. |TT1 a ferroelectric distortion along 
one of the Cartesian directions leads to a polarized struc- 
ture in the Imm2 space group, while the evolution of 
a polarization along the (111) or (111) directions leads 
to the i?3m space group. For each material we follow 
the system into its local symmetry-constrained ground 
state for each type of distortion, and we also compute 
the electric polarization in this state using the Berry- 
phase method. (Polarizations found by multiplying com- 
puted Z* values times computed displacements differ 
only slightly from those calculated using the Berry-phase 
approach.) The results are presented in the first four 
columns of Table IIIIl where the energies are reported 
relative to that of the high-symmetry FA3m structure. 
(Note that (111) denotes the direction from an A atom 
to a B neighbor, while (111) points to a B' neighbor.) 

We see that for CBTZ and KCZN the Imm2 structure 



is energetically preferred over the R3m structures, and 
the polarizations are also larger for the Imm2 structure. 
In these materials, the energies and polarizations are also 
very similar for the structures distorted along (111) and 
(111) directions. On the other hand, the PSTZ system 
behaves very differently. It is strongly polar, with the 
depth of the double- well potential, at 0.6-0.8 eV, being 
more than an order of magnitude larger than for the other 
two materials. Furthermore, the R3m structure denoted 
as (111) is now the favored structure, being significantly 
lower in energy than either the (111) R3m structure or 
the Imm2 structure. This structure also has the largest 
polarization, at 0.837 C/m^. 

We thus see that AA'BB'Og double perovskites can 
have a rich variety of polar behaviors, ranging from 
ones that remain nonpolar like KBTN, to those that are 
weakly polar like CBTZ and KCZN, and finally to the 
case of the strongly polar PSTZ. In the next subsection 
we shall see how this diversity of behaviors can be cap- 
tured in a simple analytical model. 

Before doing so, we comment briefiy on the na- 
ture of the ferroelectricity seen in these tetrahedral 
ferroelectrics. For this purpose, we have inspected 
the eigenvectors of the ferroelectric soft modes identi- 
fied in Table [TTl Letting be the sum of squares 
of soft-mode eigenvector components corresponding to 
atoms of type /i, expressed as percentages, wc find 
€ = {U,U',Cb,^b',^o) = (0.5, 45.7, 1.4, 2.1, 50.3) 
for PbSnTiZrOe, | = (0.1, 72.8, 0.2, 0.6, 26.3) for 
KCaZrNbOg, and | = (78.5, 0.5, 0.3, 0.6, 20.1) for 
CaBaTiZrOg. In all three cases, the ferroelectricity is 
found to be A-site driven, with very little involvement 
of B cations. More specifically, it is associated with a 
displacement of the smaller of the the A-site cations, 
coupled with some oxygen motion. This seems reason- 
able in retrospect, since the lattice constant of the over- 
all AA'BB'Og material will be expanded by the larger 
A atom, leaving a "rattling cage" environment for the 
smaller one. 



C. Theoretical modeling 

The results of our first-principles calculations of the en- 
ergies and polarizations of tetrahedral double perovskites 
can be modeled by expressing the energy as a polyno- 
mial in the components of the electric polarization, as 
in Landau-Devonshire theory. Symmetry considerations 
exclude certain terms in the expansion, which is then 
written as 

PiP' + Py + Pt) + ViP'Py + PyP' + P'P') + 
iP.PyP.{Pl + P^ + P!) + 0{P^) + • • • (1) 

where the energy is measured relative to that of the high- 
symmetry i^43TO space group. (Note that the 7 and ^ 
terms would vanish according to cubic symmetry.) We 
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TABLE III: Calculated polarization and total-energy reduction (relative to the high-symmetry F43m structure, per 10-atom 
cell) for distorted structures of CBTZ, KCZN, and PSTZ. Last four columns present results obtained from the models discussed 
in Sec. II V Cl (values in parentheses are exact by construction, as they were used as input to the fit). 



Ab-initio results Landau-Devonshire model 





Space 


Polarization 


Energy 


Fourth order 


Fifth order 


Material 


group 


C/m^ 


meV 


meV 


% error 


meV 


% error 


CBTZ 


i?3m(lll) 


0.137 


10.8 


11.2 


3 


(10.8) 


(0) 




i?3m(lll) 


0.136 


11.5 


11.1 


3 


11.5 







Jmm2(100) 


0.163 


16.0 


(16.0) 


(0) 


(16.0) 


(0) 


KCZN 


i?3m(lll) 


0.186 


24.9 


26.4 


6 


(24.9) 


(0) 




R3m{lll) 


0.184 


27.2 


26.0 


4 


27.5 


1 




Jmm2{100) 


0.219 


36.6 


(36.6) 


(0) 


(36.6) 


(0) 


PSTZ 


i?3m(lll) 


0.777 


681.4 


713.1 


5 


(681.4) 


(0) 




i?3m(lll) 


0.837 


794.0 


799.6 


1 


839.5 


6 




/mm2(100) 


0.709 


581.8 


(581.8) 


(0) 


(581.8) 


(0) 



denote the energy and the polarization of the orthorhom- 
bic Imm2 space group by i?(ioo) f^nd P(ioo) respectively. 
For the rhombohedral space group R3m, we have corre- 
spondingly £^(111), -P(iii), -E^{iii> and -P{iii)- -P{iii> and 
P(iii) represent the positive and negative polarizations 
of the two minima along the body diagonals. For ordi- 
nary ferroelectric materials these polarization values are 
equal in magnitude, and so are their related energies. In 
the orthorhombic case we take the z axis to be the sym- 
metry axis, so that Px = Py = and P^ = P. In the 
rhombohedral case we keep the Cartesian alignment of 
the axes, such that Px = Py = Pz = P / %/3. 

Thus, from our ab-initio calculations we have the cal- 
culated values of the six quantities -E(ioo)j 
P(iii)i ^(111): and -P(ioo) that we can use to determine 
the free parameters in Eq. ([T]). If we truncate Eq. ([T]) 
at fifth order as shown, we have six quantities to de- 
termine five parameters, thus overconstraining the solu- 
tion. Similarly, if we truncate Eq. ([1]) at fourth order, we 
overconstrain more strongly (six constraints and four pa- 
rameters). Going in the other direction to include sixth 
order in Eq. ((Ij would be problematic because several 
invariants appear at sixth order, so that more than six 
constraints would be needed to determine the system of 
equations. We therefore attempted the fits with the poly- 
nomial truncated at fourth and fifth order. When work- 
ing at fourth order, we use all values except and 
E{iii) in the fit, and then test whether we can success- 
fully predict the values of these two quantities. At fifth 
order, we omit only i and test this value from the 

fit. 

The comparison of the fitted energies with those com- 
puted from the first-principles calculations are presented 
in Table Hill Values that were included in the fit, and are 
therefore exact by construction, are shown in parenthe- 
ses. We see that even the fourth-order fit gives encourag- 
ing agreement, with a worst-case deviation of about 6%. 
However, a closer inspection reveals a specific feature of 
the fourth-order calculation that is qualitatively incor- 



rect, namely that the fitted values of i^(iii) are smaller 
than those of £:(ni) for CBTZ and KCZN, while the 
DFT calculations give the opposite trend. This relates 
to the fact that these two materials have an anomalous 
behavior in that the lower magnitude of polarization is 
associated with the deeper energy minimum when com- 
paring the and E{iii) distortions (PSTZ docs not 
show this anomalous behavior). It turns out that the 
fourth-order theory does not have enough flexibility to 
reproduce this behavior; at that level of theory, it can be 

shown that |A£;(in)| > |A£;^ni>l if l^'(iii>l > 1^(111)1, 
and vice versa. This discrepancy is removed once we go 
to the fifth-order theory; as can be seen from Table lllli 
the relative magnitudes of the energies and polarizations 
are now correct for CBTZ and KCZN, and the discrep- 
ancy between the model predictions and ab-initio calcu- 
lations improves substantially. Thus, we conclude that a 
fifth-order expansion is the minimum complexity needed 
to give a qualitatively correct description of the energy- 
polarization relations in the class of AA'BB'Oe materials 
under study here. We note, however, that there are still 
quantitative errors for PSTZ; these can only be removed 
by going to still higher order, presumably because the 
larger magnitude of polarization in PSTZ accesses higher 
terms in the Landau-Devonshire expansion. 



D. Oxygen octahedral rotations 

As we have seen in Table |TT1 most of our investigated 
materials show unstable soft modes corresponding to tilts 
and rotations of the oxygen octahedra. In our investiga- 
tion of the ferroelectric states in Sec. IIVBI we neglected 
these modes by relaxing the systems according to sym- 
metry constraints consistent with polar, but not rota- 
tional, instabilities. We now consider the effect of these 
rotations, which may compete with the ferroelectric dis- 
tortions in determining the ground state of the system. 

We use PSTZ as a test case for this purpose. As we saw 



in Table |TT1 PSTZ shows two unstable sets of modes in 
the high-symmetry F4im structure, a rotational insta- 
bility at 182icm~^ and a polar one at 160 icm^-'^. If we 
follow the polar mode distortion into the polar i?3m(lll) 
structure reported previously in Table IIIIl we find that 
it is not a local energy minimum. Instead, we find that 
this structure still has an unstable phonon of frequency 
22.7 icm~^ corresponding to rotations of the octahedra 
about the polar (111) axis (i.e., like an i?-point mode of 
the ideal 5-atom perovskite structure). Turning on these 
rotations, the system reaches a stable R3 structure at 
small rotation angles of ~0.46° and ~0.79° around (111) 
for the octahedral rotations centered on B and B' atoms 
respectively. The phonon mode that was previously un- 
stable is now found to have a positive frequency, while 
other low-lying mode frequencies and the electric polar- 
ization remain almost unchanged. This new phase is en- 
ergetically approximately equal to the i?3m(ilT) state. 
Thus, we find that while the rotations are present in the 
ground state, they are small and have little influence on 
the properties of the system. 

To check whether this structure is truly the global 
ground state, we tested what happens if we follow a dif- 
ferent route, i.e., starting from the high-symmetry F45m 
space group and following the path of the unstable octa- 
hedral rotation (mode frequency 182i in Table HI)) without 
intentionally making any polar distortion. However, be- 
cause it singles out one of the four {111} axes that were 
previously equivalent, such a rotation immediately con- 
verts the system to the polar rhombohedral space group 
i?3. Relaxation within this space group is then found 
to lead back to the same structure we found before, in 
which strong polar distortions predominate over small ro- 
tations. Thus, we again conclude that the rotations are 
unimportant for PSTZ. 

While the calculations above are specific to PSTZ, the 
symmetry analysis is more general and hints at the possi- 
bility of obtaining polar samples even if the dominant un- 
stable modes have rotational character. In fact, there is 
an interesting possibility that even if there are no unsta- 
ble polar modes in the high-symmetry FASm structure, 
an unstable rotational mode could take the system into 
a polar space group. This would correspond to the dis- 
covery of a new class of improper ferroelectrics, in which 
the primary order parameter is the antiferrodistortive ro- 
tation, but in which a polarization necessarily appears 
because the selection of a rotation axis in a parent struc- 
ture without inversion symmetry results in secondary po- 
lar distortions along that rotation axis. Moreover, it is 
interesting that electric fields could, at least in principle, 
be used to control the selection of the rotational domains 
in such a material, at least as far as selecting one of the 
rotation axes from among the four possible ones. 

We already have such cases at our disposal: both 
KSTN and SCTM have rotational soft modes at the high- 
symmetry phase, but no polar soft modes. In the case of 
KSTN, allowing rotation leads to an R3 phase which is 
about 30 meV lower in energy than the high-symmetry 
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FIG. 3: (a) Type-I and (b) type-II configuration of SCTM. 

structure, and has a polarization of 0.0023 C/m^. This 
value of polarization is indeed tiny, about 1% of that 
seen in the KCZN case, which has a comparable energy 
difference between R3m and F43m phases. In the next 
section we discuss a similar behavior that emerges from 
our studies of SCTM. 



E. Magnetic SCTM structures 

Recently, Rijnders, Blok and Blank— have succeeded 
in preparing films of the double perovskite SCTM in 
rock-salt order using layer-by-layer molecular-beam epi- 
taxy on the (111) surface of LaAlOs. To our knowledge, 
this is the first experimental realization of an AA'BB' 
perovskite system in the double rock-salt arrangement. 
While their initial characterization of this material does 
not appear to show a polar character, we were motivated 
to extend our theoretical calculations to this material in 
order to make contact with the experiments. We present 
our results on SCTM in more detail in the present sec- 
tion. 

The presence of the magnetic Mn ions makes this ma- 
terial distinct from the others studied so far. Since Sr 
and Ca are 2+ and Ti is 4+, we find Mn in its 4-1- ox- 
idation state. With its (P filling in an octahedral envi- 
ronment, this configuration is naturally compatible with 
a local spin state in which the majority t2g states are 
filled and other d states are empty. To handle the mag- 
netic nature of the Mn atom, we perform coUinear spin- 
polarized calculations, neglecting the spin-orbit coupling. 
We consider both ferromagnetic (FM) and antiferromag- 
netic (AFM) spin arrangements. The Mn atoms reside 
on an FCC lattice, which is capable of exhibiting sev- 
eral AFM structures, all of which are frustrated to some 
degree)^ We perform calculations on two of the most 
common AFM variants, the ones of type I and II illus- 
trated in Fig. [IK a) and (b) respectively. The spins are 
aligned ferromagnetically in-plane and antiferromagnet- 
ically out-of-plane with respect to (001) planes in the 
type-I structure, and with respect to (111) planes in the 
type-II structure. 

We first consider the FM spin structure, for which some 
results were also reported earlier in Tables HlHIl The 
resulting material is found to be insulating with a re- 
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TABLE IV: Lattice constants, energy (relative to FM state 
per 20-atom unit cell) and fraction of parallel nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) spins for 
three magnetic structures of SrCaTiMnOg • 



Structure 


Latt. const, 
(a.u.) 


Energy 
(moV) 


NN 


NNN 


FM 


7.15 





12/12 


6/6 


AFM (type-I) 


7.15 


-50 


4/12 


6/6 


AFM (type-II) 


7.15 


-70 


6/12 


0/6 



laxed lattice constant of 7.15 a.u. in its high-symmetry 
structure. (It also shows a rotational instability which 
we have not pursued here.) However, both of the AFM 
structures considered are lower in energy, as shown in 
Table IIVI The AFM arrangement of spins lowers the 
symmetry from F43m to tetragonal P4to2 and rhom- 
bohedral _R3to space groups for type I and II ordering, 
respectively. After relaxation, the new lattice constants 
still closely resemble that of the high-symmetry structure 
(e.g., the c/a for the pAm2 case is 0.9996). It is primarily 
the motion of the ions off the high-symmetry positions, 
brought about by the specific spin configurations, that 
defines the new space groups. Both these AFM cases are 
found to be insulators, with the type-II arrangement be- 
ing slightly lower in energy. Table ITVl also provides a con- 
cise summary of the number of parallel nearest-neighbor 
and next-nearest-neighbor spins out of all such neighbors. 

Next, we find the phonon frequencies at the F point of 
this relaxed type-II AFM cell. There is one set of three 
soft modes, but the eigenmode analysis reveals them to 
be oxygen zone-boundary phonons, ruling out the pos- 
sibility of a polar instability. This seems to be consis- 
tent with the experimental characterization of Rijnders et 
ai— However, instead of stopping here, we take the anal- 
ysis a bit further, following the spirit of the last section on 
zone-boundary rotation of oxygen octahedra. Following 
the unstable mode in thei?3m phase, we arrive at the po- 
lar i?3 space group, and find it to be lower than the R3m 
phase by 319 meV (per 20-atom unit cell). As pointed 
out before, the i?3 phase is a polar phase and exhibits an 
improper polarization induced by rotation of the oxygen 
octahedra. The calculated polarization is found to be 
very small, about 0.01 C/m^, and points along the (111) 
direction. A phonon analysis in this phase reveals that 
two modes still remain soft, indicating that the ground 
state of the system has not yet been reached. We expect 
that further relaxation along these soft modes would lead 



to a lower-symmetry structure with additional octahedral 
rotations, but still with a very small polarization. 

While it lies beyond the scope of the present investiga- 
tion, we note that it may also be of interest to consider 
the effects of epitaxial strain on the SCTM system. This 
could provide more direct contact with the experimental 
work on epitaxial growth mentioned earlier, as well as 
possibly making contact with recent work^ showing that 
SrMnOa, which is one of the parent perovskites of SCTM, 
can be driven between FM/FE and AFM/PE states via 
application of epitaxial strain. These considerations are 
left for future investigations. 

V. SUMMARY 

In summary, we have carried out a first-principles 
study of the properties of prospective AA'BB'Og per- 
ovskites having double rock-salt order. We find several 
candidate compounds that are predicted to have ferro- 
electric instabilities associated with A-site displacements, 
with PSTZ (PbSnTiZrOe) being of special interest be- 
cause of its large spontaneous polarization and peculiar 
energy landscape having four global minima along (111) 
directions and four secondary local minima along (111) 
directions. Compounds in this class may also be capable 
of exhibiting improper ferroelectricity based on rotation 
of oxygen octahedra, again stemming from the lack of 
inversion symmetry in the high-symmetry space group. 
We also predict that epitaxial strain can, under appro- 
priate conditions, induce a polarization in an otherwise 
paraelectric material of this class. The inclusion of mag- 
netic cations may provide interesting opportunities for 
the realization of novel magnetoelectric or multiferroic 
materials. 

While initial attempts at the synthesis of such ma- 
terials have not yet resulted in the demonstration of 
ferroelectricity, they do provide an existence proof that 
such synthesis is possible. We hope that the crossing of 
this hurdle will stimulate attempts at synthesis by other 
groups, and that eventual success will lead to novel ferro- 
electric materials having interesting and potentially use- 
ful physical properties. 
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